Lateral diffusive spin transport in layered structures 
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A one dimensional theory of lateral spin-polarized transport is derived from the two dimensional 
flow in the vertical cross section of a stack of ferromagnetic and paramagnetic layers. This takes into 
account the influence of the lead on the lateral current underneath, in contrast to the conventional 
ID modeling by the collinear configuration of lead/channel/lead. Our theory is convenient and 
appropriate for the current in plane configuration of an all-metallic spintronics structure as well as 
for the planar structure of a semiconductor with ferromagnetic contacts. For both systems we predict 
the optimal contact width for maximal magnetoresistance and propose an electrical measurement 
of the spin diffusion length for a wide range of materials. 
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Spintronics promises an increase in the logical express- 
ibility of electronic circuits and the integration of non- 
volatile magnetic mcmoryi^ . A quantitative and yet 
simple theory of spin transport is essential for inter- 
preting experimental results and designing practical lat- 
eral devices which are imperative in large scale integra- 
tions. The most popular analytical approach to spin 
transport extends drift-diffusion equations to include the 
spin degree of freedom^iSi&. One-dimensional equa- 
tions of this kind have successfully explained the current- 
perpendicular-to-the-plane giant magnetoresistance^ and 
the spin injection from a ferromagnetic metal into a semi- 
conductor in a vertical geometrji^iSiLSiS. Theories of spin 
transport for current in plane focus on boundary scatter- 
ing to explain the magnetoresistive effect in ultrathin lay- 
ered systemiSiii. However, many of the experiments in- 
volving non-magnetic metalsi2ii£ or semiconductorsi^ii 5 . 
are performed in the planar geometry where the thick- 
ness of the layers exceeds the mean free path. For these 
systems, application of either of the existing ID models 
is a rough approximation, unable to account for some ob- 
served features whose interpretation required a numerical 
3D treatment. The role of the contact width in length- 
ening the effective path between two leads versus the spin 
diffusion length was mentioned^. The effect of leakage 
from two-dimensional electron gas into the ferromagnetic 
gate of finite lateral size on the spin currents in the chan- 
nel has been investigated^. What remains lacking is a 
simple theory in which the spin injection from contacts 
of finite width is normal to the spin current in the plane 
of the conduction channel. 

In this paper, starting from a study of the two dimen- 
sional current flow in the vertical cross section of the 
planar system, we derive an analytical ID theory of spin 
transport in lateral structures. The analysis takes into 
account the lateral extent of the interface between ferro- 
magnetic (FM) and normal (N) materials. We obtain a 
transparent physical picture and derive useful analytical 
formulas. The most important result obtained is the de- 
pendence of the magneto-resistance (MR) on the contact 
width. The knowledge of this aspect of lateral spin- valve 
physics is essential to the design of structures with high 
MR. In addition, our analysis suggests a technique to ob- 
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FIG. 1: A two-dimensional sketch of a planar system with 
ferromagnetic and normal materials. A-D represent bound- 
aries which may be connected to the external circuit. All of 
the labeled lengths enter into our ID transport equations. 



tain the spin diffusion length of a paramagnetic material 
by a simple set of electrical measurements. For com- 
pounds, such as silicon, which are inaccessible to optical 
characterization of spin polarizationi^^, development of 
the electrical method is of great importance. 

Fig. depicts a typical planar system, homogeneous 
in the z direction. A model of one dimensional current 
flow from lead B via the normal region to lead C cannot 
describe the flow under the leads because of the simulta- 
neous in-plane flow and the injection or extraction from 
above. Hence, we must first consider the two-dimensional 
transport governed by the spin diffusion equation: 



V 2 /i s (x,y) = 



Hs{x,y) ~ Hs{x,y) 



(1) 



where fi s is the spin dependent electrochemical potential 
with s=± denoting the spin species. The characteris- 
tic distance of spin flip is L s = ^/D s r Sj _ s , where r Si _ s 
is the spin-flip time and D s is the diffusion constant of 
the spin component s. The measurable spin diffusion 
length L, on whose scale the splitting of fi changes, is 
given by L~ 2 =L~ 2 + LZ 2 . These equations hold for fer- 
romagnetic and normal metals 5 and also hold for non- 
degenerate semiconductors with the additional condition 
of a small electric field^. We neglect interfacial spin scat- 
tering and use Ohm's law across the interfaces. Since the 
current is driven by the gradient of the electrochemical 
potential, the continuity of the normal component of the 
spin current across two adjacent layers i and k provides 
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the boundary conditions, 
at (V • VmJ) = Gf( Ms fc ~ Ml) = -°s (n k • V M *) , (2) 

where ct* and n 1 refer to layer i conductivity and out- 
ward interface normal, respectively, and G 1 ^ is the spin- 
dependent barrier conductance. Although Fig. ^ illus- 
trates only two layers (normal and FM), the equations 
derived below hold for multilayers. 

To reduce the essential features of the two dimensional 
flow to one dimension, we define the vertical (y) average 
of fj, s in each layer, 



1 f y i 



dy^ s (x,y), 



(3) 



with h l being the thickness of the layer between its 
boundaries, (yQ,y[). For thin layers, fj, % s (x,y) may be 
replaced by its vertical average (x) . The conditions of 
validity follow the requirement that the gradient correc- 
tion along y to fi s in Eq. @ be negligible to 0(h 2 ), 



(i) h -C <j s /G s , and (ii) h < L s 



(4) 



with the help of Eqs. and J3J). Under these conditions 
we can derive a set of ID equations governing the lateral 
(i.e., in the layer plane) spin transport. Integrating out 
the y dependence in Eq. and using Eq. |3J) yields, 



t+1 r 



dx 2 



(LI) 2 



aih 1 



C(x)-C +i (x) 



(5) 



This kinetic theory assumes that the relevant length 
scales exceed the electron mean free path. Thus, the 
present analysis cannot address the quantum realm of 
the MR effect in ultrathin layered systems^. 

We now divide the lateral transport into regions of 
vertical stacks. Consider a region with a vertical stack 
of Nl layers. For example, in the region of width wl 
covered by the left lead in Fig. ^ Nl — 2, excluding the 
insulator. Transport in the layers of the stack is governed 
by 2Nl coupled members of Eq. (JSJ which in the matrix 
form are simply, 



(d 2 /dx 2 )$, = M£. 



(6) 



with the column vector £ of elements Q and the positive 
definite matrix M of elements from Eq. J5J) . The solution 
has the form, 



£(x) = l(a + bx) 



Nl-1 

E • 

3=1 



j [Pi* 



(7) 



of the system and to each other is maintained by the 
boundary conditions. Between stacks they are given by 
the continuity of £ s (x) and their first derivatives (cur- 
rents) through each homogeneous layer. At the outer- 
most boundaries (like A and D in Fig. ^) the conditions 
are prescribed by the external driving terms of the elec- 
trochemical potential in the form of either a constant 
voltage maintained at an electrode or an injection cur- 
rent (including zero) through an interface. Applying cur- 
rents/voltages to B and C interfaces results in the pres- 
ence of inhomogenous terms in Eq. 10. These boundary 
conditions provide a unique set of solutions for the pa- 
rameters a,b,pj,qj in Eq. J7J). This method of solution 
is a considerable simplification compared to solving the 
2D spin diffusion equation for the same planar structure, 
and yet it includes the influence of the vertical layers on 
the lateral spin current flow. In the following, we will 
show the versatility of this simplified approach as well as 
testing its validity against the 2D solution. 

The first illustrative example involves a semiconductor 
in the normal layer of the structure in Fig. ^ A and D 
are disconnected from outside while B and C are voltage 
biased. The spin diffusion length in the normal layer is 
Ljv. We divide the normal channel into sections belong- 
ing to 5 vertical stacks as shown in Fig. ^ In the middle 
section of the channel, the solution is£*&, 

£±(a0 = (e/a N )Jx + fi±( Pc e x / L " +q c e- x / L *), (8) 

where J, [i,p c and q c are constants to be determined by 
the boundary conditions. The total current flowing be- 
tween the two leads is / = Jh N per unit length of the 
structure in the z direction. For the two sections outside 
the footprint defined by the two leads, we introduce the 
notion of the "open" versus "confined" geometry depend- 
ing on whether cIl , dn 3> or <C ijv • I n both geometries 
the pattern of the net charge current occurs only between 
the leads B and C. This means that the outside sections 
lack the linear term of Eq. JSJ. Additionally the "open" 
geometry includes only an exponentially decaying solu- 
tion away from the closer lead, and in contrast with the 
confined geometry, spin polarization extends noticeably 
outside of the footprint between B and C, thus reducing 
the spin accumulation densities. In each of the vertical 
stacks containing a lead, due to vastly different conduc- 
tivities of the semiconductor and ferromagnetic metal we 
can decouple the equations for £ s in the N semiconductor 
and in the FM metal and solve for the eigenvalues of M 
in the former, 



[a + l± ^T^]/(2L 2 N ), 
(a, (3) = 2L 2 N (G + ±G^/(a N h N ), 



\ 2 



(9) 



where 1, a column of ones, is an eigenvector of M with 
zero eigenvalue and vj is an eigenvector with eigenvalue 
X 2 . The connection of the vertical stacks to the outside 



where the first of each pair of symbols (s, c) or (a,/3) 
takes the upper sign. The electrochemical potential in 
the FM is practically a constant (£, FM =—eV), and in 
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FIG. 2: (a) Magneto-resistance effect versus the contact width 
of a GaAs channel at 300K. The solid (dashed) line denotes a 
confined (open) structure. The dots are the results of a 2D nu- 
merical computation in the confined geometry, (b) Magneto- 
resistive effect versus a for three cases of spin selectivity, j3/a 
in the confined geometry, calculated for the optimal contact 
width corresponding to the same a value shown in (c). 



the semiconductor layer under the biased lead is 
£ ± (x) = -eV+(l±X) 
+ (ATl 



pe + qe 



re x " x + se- x ' x 



(10) 



where p, q, r and s are constants to be determined by 
boundary conditions and A = cot[i tan -1 0\. Consider 
the case that spin polarization is robust so that a and 
(3 are comparable. If a <C 1, then A c <C and 
X s ~ l/Ljy. The s-mode is limited by the spin diffu- 
sion constant, and it corresponds to spin accumulation 
(A»l in this case). If a 3> 1, then both eigenvalues are 
nearly independent of Ln, and neither of the eigenvectors 
is a pure spin mode A~l: the inhomogeneity of injection 
dominates the spatial dependence. 

To illustrate the effect of the contact width on MR, 
we use the experimentally verified barrier parameters of 
a Fe/GaAs system at 300Ki2i3I. The tunneling conduc- 
tances are of the order of 10 2 -10 3 ri~ 1 cm -2 with the ratio 
of spin-up to spin-down conductance G+/G_ ~ 2. The 
N channel doping is no=4 x 10 15 cm~ 3 and the spin relax- 
ation time is t s —80 ps<2i which corresponds to L N ~ 1 /an. 
The thickness of the N channel is h N — 100 nm, and the 
inner channel length is <i m =200 nm. 

The MR effect is defined as (I P -I AP )/I P with I p 
(I AP ) being the total current when the FM contacts 
are magnetized in parallel (antiparallel) directions. The 
calculated MR is shown in Fig. Ufa) as a function of 
the contact width, where for simplicity we have used 
w=Wl=Wr. In the confined geometry, dL=dji=0, the 
effective ID method (solid line) is compared with the 
results of 2D numerical calculation (dots), showing ex- 
cellent agreement. We note the deleterious effect of the 
open geometry (dashed line) compared with the confined 
geometry (solid line) because of the weaker spin accumu- 
lation in the semiconductor channel which produces the 



MR effect. 

Note the existence of the optimal contact width for 
the maximum MR effect which arises out of the balance 
between spin injection and spin relaxation in the chan- 
nel. The spin accumulation in the normal conduction 
channel is built up from injection along the width of the 
ferromagnetic contact. When the width is small, inad- 
equate build-up results in small spin injection and spin 
accumulation as a percentage of the equilibrium carrier 
density. Consequently the MR effect is small. However, 
when the contact width is too large, its resistance be- 
comes very small, the MR effect in the N channel is again 
negligible as in the conventional one-dimensional model—. 
Finally, when the contact width w exceeds the spin dif- 
fusion length Ln, the build-up from vertical injection 
along the width of the contact beyond Ln becomes in- 
effective, the spin injection and the MR effect approach 
asymptotic values. The difference between the open and 
confined structures is removed, as can be seen from the 
merging of the solid and dashed lines in Fig. Ufa) . The 
magnitude of MR at the optimal contact width depends 
also on the spin selectivity of the barriers, f3 in Eq. 
as illustrated in Fig.|2Jb). Fig.^c) shows a much weaker 
dependence on the spin selectivity of the optimal contact 
width in units of the spin diffusion length (w op t/L N ). 
The ratio depends only on a as defined in Eq. 10. For 
ci»l, w opt /L N ~ 2/ a. For a<l, the optimal contact 
width w opt /L N is approximately 6/5 or 3/8 for the open 
or confined geometry, respectively. These results enable 
the extraction of the spin diffusion length of a test ma- 
terial by measuring the MR of several structures with 
different contact widths from the same growth. 

To demonstrate the use of our method for the all- 
metallic system, we turn to a second example. In the 
planar structure of Fig. ^ a constant current is injected 
from region A into a 10 nm paramagnetic metal layer. 
The structural parameters in terms of the spin diffu- 
sion length of the N layer are 0. 2e?£=0. 2du=2dj\i=LN- 
The FM layers B and C are of the same thickness 
and serve as floating contacts. We use the Al/Co 
system parameters at 300K where ratio of spin-up to 
spin-down conductance is G+/G_=2 23 and their sum is 
1.5 x 10 9 £l -1 cm~ 2 . The normal layer intrinsic parame- 
ters are <7 N =3x 10 5 f2 _1 cm _1 , L N =0.35 /im. In the Co 
layers the majority channel conductivity and spin diffu- 
sion lengths are, respectively, O.ler^ and 12 nm. The 
respective values of the Co minority channel are 0.04<r iv 
and 20 nm. For ease of toggling between P and AP con- 
figurations we use a width ratio of 2 between left and 
right Co regions. For the FM/N regions the solutions are 
given by Eq. @. In the middle N channel the solution 
of £ s (x) is given by Eq. © while at the outer channels 
the form of £ s (a;) lacks the diverging exponential. 

For each magnetic configuration (P or AP) we calcu- 
late the voltage drop Vbc with the electrochemical poten- 
tials averaged along each Co layer. The MR is defined as 

li-^BcABcfl- Fi g- H a ) shows the MR effect - lt P eaks 
at the right contact width of wr ~ 0.5L N . The opti- 
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FIG. 3: (a) The magneto-resistive effect versus the right con- 
tact width in Al/Co system at 300K. (b) The current density 
profile in the Al layer for parallel magnetic configuration at 
the optimal contact width. Solid (dash) lines denote spin-up 
(down) current components. The total current is constant 
where in the dip regions along the Al/Co interfaces currents 
propagate in both layers, (c) Same as (b) for the anti-parallel 
magnetic configuration. 

mal contact width mostly on a, insensitive to varying 
the spin selectivity, the middle N channel length, or the 
outer N channel lengths (not shown) . Fig. |3fb) and (c) 
show, respectively, the current densities in the normal 
layer for parallel and antiparallel magnetic configuration 
at the optimal contact width. The solid (dashed) lines 
depict spin up and spin down current components along 



the x direction. Despite the small MR effect, its observa- 
tion is well within existing experimental abilities. Thus a 
similar analysis to the previous FM/SC example leads to 
an electrical method to measure the spin diffusion length 
of the normal material. However, the extraction of Ln 
requires that it greatly exceeds the FM spin diffusion 
length. 

In summary, we have presented an effective ID theory 
which describes spin dependent transport beneath and 
between ferromagnetic contacts in the lateral geometry. 
Under the conditions of thin layers and non-ohmic inter- 
faces, the derived set of coupled linear equations govern- 
ing the lateral diffusive spin currents is almost as accurate 
as the 2D spin diffusion equation but much simpler to use. 
This method retains the important role of the contact 
width, predicting an optimal contact width for the maxi- 
mal magneto-resistive effect in the lateral spin- valve. We 
have also quantified the effect of "open" channels, where 
spins can diffuse away from the active region underneath 
the FM contacts, in contrast with the "confined" struc- 
tures, where the spin accumulation is kept between the 
two contacts. Our results provide crucial guidelines for 
design of lateral spin-valve type devices. They also yield 
a method of extracting the spin-diffusion length in an 
all-electrical measurement. At room temperature where 
the spin diffusion length lies in the sub-micron region, 
this method does not have the difficulty of wavelength 
resolution in the optical techniques. 
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